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1.  Introduction 


Composite  materials  are  currently  used  for  lightweight  pressure  vessels  and  highly  efficient 
rotors  for  energy  storage.  For  both  applications,  pre-stresses  are  built  in  during  fabrication  of  the 
cylinders  through  a  “press-hf  ’  procedure  to  enhance  the  mechanical  performance.  For  pressure 
vessels,  this  pre-stress  is  achieved  by  pressurizing  the  vessel’s  interior  liner  prior  to  winding  the 
composite  overwrap.  Therefore,  after  manufacturing,  the  liner  is  in  a  state  of  hoop  compression 
and  the  composite  overwrap  is  in  tension.  For  the  rotor  application,  the  rotors  are  subjected  to  a 
radial  compression  prior  to  operation.  Accordingly,  the  centrifugal  force  resulting  from  the 
rotation  of  the  rotors  generates  tensile  stresses  in  the  radial  and  circumferential  directions.  Since 
the  composite  rotors  are  mainly  circumferentially  reinforced  (filament-wound  cylinders),  the 
radial  tensile  stress  is  critical  to  the  ultimate  perfonnance  of  the  rotors.  It  is  essential,  then,  to 
design  and  build  the  rotors  with  radial  precompression.  However,  polymer  matrix  composites 
generally  creep  over  a  long  period  of  time,  especially  at  an  elevated  temperature  (1 ).  The 
associated  stress  relaxation  in  the  composite  will  result  in  the  loss  of  the  pre-stresses  and  lead  to 
a  potential  failure.  The  objective  of  this  investigation  is  to  develop  an  analytical  method  to  study 
the  viscoelastic  behavior  of  thick-walled  composite  cylinders.  The  analysis  can  be  applied  to  the 
design  of  flywheel  machinery  and  composite  pressure  vessels. 

To  date,  activities  in  the  research  of  viscoelasticity  have  mainly  concerned  isotropic  materials, 
including  studies  by  Muki  and  Sternberg  (2),  Schapery  (3),  Williams  ( 4 ),  and  Christensen  (5). 
These  basic  theories  of  viscoelasticity  were  then  extended  to  the  area  of  heterogeneous  and 
anisotropic  materials  for  a  variety  of  applications.  Hashin  (6)  used  the  effective  relaxation 
moduli  and  creep  compliances  to  define  the  macroscopic  viscoelastic  behavior  of  linear 
viscoelastic  heterogeneous  media  and  its  implementation  in  viscoelastic  modeling.  Schapery  (7) 
examined  the  applications  of  the  correspondence  principle  and  the  general  formulation  of  linear 
viscoelastic  boundary  value  problems  of  composite  materials,  including  the  thennal  viscoelastic 
problems  for  thermorheologically  simple  materials..  Rogers  and  Lee  (3)  investigated  the 
viscoelastic  behavior  of  an  isotropic  cylinder.  Tzeng  and  Chien  (9)  used  a  multiscale  modeling 
modeling  technique  to  study  heterogenous  viscoelastic  solids  containing  evolving  cracks.  In 
addition,  finite  element  packages,  such  as  ABAQUS,  ANSYS,  and  DYNA3D,  are  not  suitable 
for  the  viscoelastic  analysis  of  composite  cylinders  because  of  the  lack  of  anisotropic  viscoelastic 
elements.  A  viscoelastic  analysis  was  proposed  by  Chien  and  Tzeng  (10)  to  solve  thennal  stress 
relaxation  in  an  anisotropic-laminated  composite  cylinder.  Souza  and  Allen  (11)  used  a 
multiscaled  modeling  technique  for  solving  viscoelastic  solid  with  microcracks.  An  extensive 
review  on  viscoelastic  solids  by  Wineman  (12)  illustrated  the  recent  progress  in  this  particular 
field. 


1 


2.  Viscoelastic  Formulation 


The  following  research  studied  the  quasistatic  viscoelastic  behavior  of  the  liner  in  a  thick- 
laminated  composite  cylinder  with  thermal  stress  caused  by  an  elevated  temperature  change. 
Accordingly,  the  field  temperature  is  constant.  The  analysis  accounts  for  ply-by-ply  variation  of 
properties,  temperature  changes,  and  fiber  orientations.  The  thick  cylinder  is  assumed  to  be  in 
the  absence  of  thermomechanical  coupling  and  in  a  state  of  generalized  plane  strain  where  all 
stress  and  strain  components  are  independent  of  the  axial  coordinate  (9,  10).  Moreover,  due  to 
the  nature  of  axisymmetry,  all  stress  and  strain  components  are  also  independent  of  the 
circumferential  coordinate.  The  mechanical  responses  of  this  thick  composite  cylinder  will, 
therefore,  have  to  satisfy  only  the  governing  equation  in  the  radial  direction.  Invoking  the 
Boltzmann  superposition  integral  for  the  complete  spectrum  of  increments  of  anisotropic 
material  constants  with  respect  to  time,  the  thermoviscoelastic  constitutive  relations  of  the 
anisotropic  composite  cylinder  can  be  derived  in  integral  fonns.  Since  the  thick  composite 
cylinder  is  subjected  to  a  constant  elevated  temperature  and  boundary  conditions  are  all 
independent  of  time,  formulations  of  the  linear  thermal  viscoelastic  problem  can  have  forms 
identical  to  those  of  the  corresponding  linear  thermoelastic  problem  by  taking  advantage  of  the 
elastic-viscoelastic  correspondence  principle.  In  other  words,  all  of  these  integral  constitutive 
equations  reduce  to  the  algebraic  relations,  which  are  very  similar  to  those  developed  for 
thennoelastic  media  when  they  are  Laplace-transformed  by  means  of  the  rule  for  convolution 
integrals.  The  thennoelastic  analysis  can  thus  be  used  to  derive  the  transformed  thermal 
viscoelastic  solutions  in  the  frequency  domain. 

The  Boltzmann  superposition  integral  of  stress,  cr,/  (i,j  =1,2,  3),  and  strain,  sy  (i,j  =1,2,  3), 
relation  for  an  isothermal  viscoelastic  problem  with  a  constant  temperature  increases  AT,  and  the 
thennal  expansion  coefficient  o,u(  T)  is 


cJij(t)  =  \'o  Cy(T,t  -  t)  pj(T,t)AT , 

where  C f  (T,  t)  is  the  relaxation  modulus  dependent  on  temperature  T  and  time  t, 


Cy  (T(t),t)  =  C*J(T0,W) 


and 


Mt)  = 


t 

ClT  (T(t))  ' 


(1) 

(2) 

(3) 
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Here,  To  is  the  base  temperature  and  ar  is  the  temperature  shift  factor.  fb/Tj)  is  given  by 
/3ij(T,t)  =  Cijkl(T,t)»a.ki(T).  It  is  often  desirable  to  use  the  inverse  form  of  the  constitutive  relation 
in  equation  1 , 


,(t)  =  fo  A( T’t  -  r) 


d<?ki(  r) 
dr 


dr  +  y/JT,t)AT, 


(4) 


where  WjfT.t)  is  the  tensor  product  of  the  creep  compliance  Ay  (T,t)  and  the  thermal  creep 
coefficient  (pkl .  Since  the  elevated  temperature  change  AT  is  constant  above  some  reference 
value  in  time,  the  relaxation  moduli  and  creep  compliance  are  evaluated  at  that  reference 
temperature,  regardless  of  whether  or  not  the  material  is  thennorheologically  simple,  by 
employing  the  temperature  shift-factor. 

The  Laplace  transfonn  of  a  function /(f)  is  defined  as 

7  =  f(s)  =  C  e"f(t)dt ,  (5) 

where  s  is  the  Laplace  transform  variable.  Accordingly,  the  temperature  is  considered  as  a 
constant.  Applying  equation  3  with  the  convolution  rule  to  equations  1  and  2  reduces  the 
integral  constitutive  equations  to  the  following  algebraic  relations: 


and  the  inverse  form 


~  ki  ~  AT 

=  C  i  '  fiv~ 


respectively,  where 


kl  — 

ij  <7] 


+  ¥ 
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C  y  =  S  dy  ,  A 


=  s  Ay ,  a  =  Say $  ,=  s#  fi  kl=  C  % 
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and 


Furthermore,  it  can  be  shown  that 


~  ~A  a  r 
Y  kl=  A  0  •  <f>  y  . 


A  k! 

— 

c  k! 

v  _ 

V  _ 
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(8) 
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(10) 


Consider  a  filament-wound,  axisymmetric,  thick  composite  cylinder  consisting  of  N  layers  with 
the  axial  coordinate  z,  the  radial  coordinate  r,  and  the  circumferential  coordinate  9,  as  shown  in 
figure  1.  The  composite  cylinder  has  the  inner  radius  a,  the  outer  radius  b,  and  the  length  L. 
There  is,  therefore,  a  corresponding  thermoelastic  problem  with  the  transformed  displacement 
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Figure  1.  Cylindrical  coordinate  system  (r,  9,  z)  is  used  to  derive  the  governing  equations 
for  the  laminated  composite  cylinder.  Radial  displacement  and  stress 
continuity  conditions  are  imposed  to  obtain  the  multilayered  cylinder  solution. 


components  u,  v,  and  w  in  the  axial  direction,  the  circumferential  direction,  and  the  radial 
direction,  respectively,  in  each  layer.  The  axisymmetric  character  of  the  thick  composite 
cylinder  and  the  assumption  of  the  state  of  generalized  plane  strain  leads  to  a  simplified 
displacement  field,  which  reflects  the  circumferential  independence  and  only  radial  dependence 

of  w : 

u(r,0,z)=u(r,z), 
v(r,  0,  z)  =v(r,  z), 
and 

w(r,  0,  z)  =w(r). 

Since  each  layer  of  the  thick-laminated  cylinder  is  cylindrically  monoclinic  with  respect  to  the 
global  coordinates,  there  is  no  coupling  between  transverse  shears  and  other  deformations.  The 
vanishing  shear  traction  boundary  conditions  and  interface  continuity  conditions  thus  generate 
zero  out-of-plane  shear  traction  and  shear  strains  for  each  layer.  Moreover,  owing  to  the  absence 
of  torsional  deformation,  the  transformed  displacement  components  u  and  v  become 


(ID 

(12) 

(13) 
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—o 
£  Z 


(14) 
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and 


(15) 


where  the  constant  quantity  s  has  the  physical  interpretation  of  transformed  axial  strain  of  a 
layer.  In  fact,  s  ,  according  to  the  present  formulation,  also  represents  the  transformed  axial 
strain  of  the  entire  composite  cylinder.  The  calculation  of  g  requires  the  knowledge  of  end 
boundary  conditions  and  will  be  given  later.  Likewise,  solving  for  w(r)  requires  the  information 
of  transformed  strain  components,  the  constitutive  equations,  as  well  as  the  equilibrium 
equations. 


The  previously  transfonned  displacement  field  gives  the  transformed  strain  components  in 
cylindrical  coordinates: 


d  w(r)  - 
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The  unabridged  form  of  the  constitutive  equation  (4)  for  each  layer  in  cylindrical  coordinates 
with  the  radial  coordinate  r  normal  to  the  plane  of  symmetry  is  expressed  as 
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Furthermore,  from  the  previous  discussions,  it  can  be  shown  that  two  of  the  three  equilibrium 
equations  are  satisfied  automatically.  The  only  nontrivial  equilibrium  equation  is  the  one  in  the 
radial  direction: 
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d  (7 IT  _|_  ( y n 

dr 


CT  98 


=  0. 


(21) 


Substituting  equations  6,  7,  and  8  into  9,  the  transformed  stress  components  cr „  and  <j0o  are 
obtained  in  terms  of  the  transformed  radial  displacement  w .  Incorporating  the  resulting 
cjrr  and  (jun  functions  with  equation  10  gives  a  nonhomogeneous  Euler  differential  equation  of 

w  for  a  layer 


,  d2w  dw  r 

r — y  +  r—  -  Aw  =  — 
dr  dr  C 


where 


33  L 


AT  _  _  ^ 

(  B rr  "  Hee)  (  Cl3~  C 12)  £ 


(22) 


Solving  (1 1)  for  w  yields 


rj  _  c„ 

A  c,/ 

(23) 

=  +  A2rA-  +  wp. 

(24) 

wP  =  f~s°r  +  f3r. 

(25) 

J  _  Cn  *  C/i 
’  C33  -  C22  ’ 

(26) 

S 

f  3  = 

3  C33  -  C22 

(27) 

A  2V  ~  l 

S  -  ^7  [  B,r  -  B////J  ’ 

(28) 

and  A/  and  A2  are  coefficients  to  be  determined  from  boundary  and  continuity  conditions. 
Finally,  it  is  understood  that  the  initial  condition  of  the  original  thermoviscoelastic  problem  is 
displacement-free  state  of  rest.  The  boundary  condition  is  of  free  traction  and,  hence,  of  free 
transfonned  traction  on  both  inner  and  outer  circular  surfaces: 


crrr  ~  <7g,.  ~  <7,r  ~  0  at  r  -  a,  b . 

On  both  end  surfaces,  stress  resultants  are  zero: 

n  r 


Z  J  <T=rdr  =  =  oz0  =  0  at  z  =  0,  L. 

k= 1  n 


(29) 


(30) 
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The  r,-  and  rQ  are  inner  and  outer  radii,  respectively,  of  the  Mi  layer.  The  continuity  conditions  at 
each  interface  between  two  adjacent  layers  require  continuous  radial  traction  and  continuous 
radial  displacement  at  any  instant  as  shown  in  figure  1 .  Thus,  when  written  in  the  transformed 
form,  they  become 

~(k) 

CJ rr,o 

and 

w, 

where  k  =  -  1;  subscripts  i  and  o  denote  inner  and  outer  surfaces,  respectively. 

Accordingly,  the  formulation  accounts  for  ply-by-ply  variations  of  material  properties  and 
temperature  change.  The  matrix  fonn  numerical  solution  procedure  with  parallel  computing 
techniques  resolved  the  complexity  and  time-consuming  calculation  procedures  in  Laplace 
transfonn  of  a  multilayered  composite  cylinder  (10). 


—(k+l)  _  n 
(7n,i  ~  0 


(k+l) 

=  W.- 
VV  >1  ? 


(31) 

(32) 


3.  Experimental  Validation 


To  validate  the  viscoelastic  model  described  above,  a  test  article  was  fabricated.  This  test  article 
consisted  of  a  steel  cylinder  liner  overwrapped  with  a  carbon-fiber  epoxy  shell  using  a  high- 
tension,  wet-filament  winding  process.  The  inner  diameter  (ID)  of  the  steel  cylinder  was  4.72  in 
(120  mm),  the  wall  thickness  12.50  mm,  and  the  length  1  m.  The  fiber  used  was  Hexcel  IM9,* 
12,000  filaments  per  tow,  and  the  epoxy  system  was  Huntsman  EPON  862  with  an  amine 
curative  (bis(p-aminocyclohexylmethane),  PACM). 

A  2-m  length  of  IM9  carbon  fiber  was  weighed  ahead  of  time  to  detennine  its  cross-sectional 
thickness  by  dividing  the  mass  per  length  by  the  material  density  (manufacturer’s  data  sheet 
reports  a  density  of  1.80  g/cm  ).  Using  this  approach,  the  cross-sectional  area  of  a  dry  IM9  tow 
was  calculated  to  be  0.0003  in"  (0.187  mm  ).  This  cross-sectional  area  was  used  to  determine 
appropriate  filament-winding  parameters  to  process  the  test  article  with  the  desired  architecture, 
thickness,  and  length.  Three  tows  of  IM9  carbon  fiber  were  wet-filament- wound  with  a 
programmed  bandwidth  of  0.15  in  (3.81  mm)  and  a  fiber  architecture  of  [±842/902]s.  Note  that 
each  ±84°  ply  is  essentially  comprised  of  one  +84°  ply  and  one  -84°  ply  that  are  helically 
“meshed”  together  as  they  are  deposited  during  winding.  As  such,  the  designed  layup  is  actually 
30  plies  thick.  This  architecture  was  chosen  with  the  assumption  of  a  65%  fiber  volume  fraction 
in  order  to  achieve  a  final  wall  thickness  in  the  composite  of  0.54  in  (13.65  mm).  The  ±84° 


Hexcel  is  a  registered  trademark  of  Hexcel  Corporation,  Stamford,  CT. 
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helical  plies  were  incorporated  to  stabilize  the  wet  composite  material  on  the  mandrel,  thus 
preventing  the  axial-direction  “extrusion”  of  ungelled  underlying  material  that  otherwise  occurs 
(especially  under  high-tension  winding)  on  pure  90°-thick  architectures. 

The  70-cm-long  carbon  fiber  overwrap  was  filament-wound  onto  the  1-m-long  steel  tube  using  a 
McClean  Anderson  4-Axis  SuperHornet  filament  winder.  The  winder  was  adapted  with  a  high- 
tension  apparatus  (designed  and  built  at  the  U.S.  Army  Research  Laboratory  [ARL],  Aberdeen 
Proving  Ground,  MD)  to  fabricate  the  part  with  an  intentionally  elevated  level  of  residual  stress. 
The  tension  is  electronically  controllable  and  was  set  at  a  value  of  85  lb  (387  N)  for  the  entirety 
of  the  part.  Assuming  a  fiber  volume  fraction  of  65%,  this  value  of  tension  corresponds  to  a 
fiber  direction  stress  level  of  65  ksi  (450  MPa).  Typical  filament  winding  is  performed  with  a 
nominal  level  of  fiber  direction  stress  (say,  3-4  ksi) — just  enough  to  keep  the  fibers  straight  and 
to  consolidate  the  layers. 

After  winding,  the  part  was  cured  at  121  °C  for  8  h  in  a  convection  oven.  Following  the  cool¬ 
down,  the  part  was  cut  into  22  ring-shaped  specimens,  each  1  in  (2.5  cm)  in  axial  length,  as 
shown  in  figure  2.  Two  specimens  were  instrumented  with  hoop-oriented  strain  gages  placed  on 
the  steel  ID  and  tested  immediately  (with  no  postcure  thennal  conditioning)  by  cutting  the 
composite  sheath  off  the  steel  ring  and  measuring  the  strain  growth.  Four  additional  specimens 
were  thermally  conditioned  by  soaking  for  10  days,  two  specimens  at  95  °C  and  two  at  1 15  °C. 
After  the  thermally-soaked  specimens  cooled,  strain  gages  were  installed  on  the  ID  of  the  steel 
and  sheaths  were  cut  off.  Measured  strains  for  all  specimens  are  summarized  in  figure  3.  The 
measured  strain  can  be  used  to  calculate  the  compressive  stress  relieved  that  causes  the  ring 
expansion.  Furthermore,  this  compressive  stress  is  equal  to  the  radial  compressive  stress  at  the 
interface  between  the  composite  and  steel  prior  to  the  sheath’s  cut-off.  Accordingly,  the 
viscoelastic  analysis  is  used  to  calculate  the  stress  in  the  overwrapped  cylinder  subjected  to  the 
temperature  condition.  The  developed  viscoelastic  analysis  is  used  to  calculate  stress  profile  and 
compared  with  the  compressive  stress  resulting  from  the  ring  test  proposed  in  figures  2  and  3. 


Figure  2.  Illustration  of  parent  tube  and  ring  specimen. 
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Figure  3.  Results  from  ring  relaxation  experiment. 

The  viscoelastic  properties  of  the  composite  matrix  (EPON  862)  were  characterized  using 
dynamic  mechanical  analysis  (DMA)  (Q800,  TA  Instruments,  New  Castle,  DE).  Beam  samples, 
nominally  2  x  0.5  x  0.125  in  (50  x  12.5  x  3  mm),  were  tested  in  a  DMA  fitted  with  a  1.37-in 
(35-mm)  dual-cantilever  beam  fixture.  Experiments  at  a  range  of  temperatures  25-1 15  °C  were 
performed  in  stress  relaxation  mode  for  20  min  and  -0.2%  strain.  The  master  curve  of  relaxation 
modulus  at  various  temperatures  and  the  time  temperature  shift  factor  are  illustrated  in  figure  4. 
The  master  curve  can  be  curve-fitted  into  a  function  to  represent  the  effect  of  time  and 
temperature  of  composite  matrix.  The  transverse  and  shear  compliances  of  composite  at  25  °C 
(room  temperature)  are  illustrated  in  the  following  forms: 

s2  (2  t=)  s°2  («)ao  , 

s6(6t=)s°6  OO00  , 

S22  =  7.5328x1  O'7/  psi, 
and 

S°66  =  1.3834  xlO'6/  psi 

The  compliance  in  the  fiber  direction,  Si  =j5  .  x  p  o/f>  ,  A)’.?(0=  and  Poisson’s  ratios, 

Vl2  =  Vl3  =  0.3,  y2J  =  0.36 ,  are  assumed  to  be  time-independent  for  this  example  simulation. 

The  analysis  of  the  ring  test  is  performed  using  the  room  temperature  properties,  as  illustrated. 
Accordingly,  time-temperature  superposition  methodology  (13)  is  used  to  adjust  the  time 
duration  for  creep  and  stress  relaxation  computation.  The  data  from  the  rings  conditioned  at 
elevated  temperatures  of  95  °C  is  used  for  comparison  with  the  model  simulation  since  115  °C  is 


(33) 

(34) 

(35) 

(36) 
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Figure  4.  Master  curves  of  creep  compliance  and  time -temperature  superposition. 

actually  above  the  glass  transition  temperature  of  the  matrix.  Based  on  the  principle  of  time- 
temperature  superposition,  the  viscoelastic  response  in  the  ring  conditioned  at  95  °C  for  10  days 
(14,400  min)  will  be  equivalent  to  a  ring  at  the  room  temperature  for  almost  infinite  time  (1010 
min),  as  shown  in  figure  4. 

Figures  5  and  6  show,  respectively,  the  radial  and  circumferential  stress  profiles  from  the 
simulation.  The  initial  stress  is  preset  to  be  equal  to  the  stress  state  in  the  composite  overwrap 
cylinder  used  in  the  experimental  validation.  The  interface  of  composite  and  steel  region  is  at 
the  radius  of  about  2.75  in,  where  a  discontinuity  of  stress  can  be  observed  clearly.  This 
discontinuity  is  attributed  to  the  difference  of  material  properties  of  composite  and  steel  liner. 
Stress  profiles  at  three  time  intervals  (initial,  two  years,  and  infinite)  are  illustrated.  The  radial 
stresses  calculated  from  the  ring  measurement  are  superposed  on  the  curves.  A  reasonable 
agreement  is  achieved.  Stress  relaxation  due  to  viscoelastic  response  of  composite  is  properly 
predicted  as  a  function  of  time.  The  hoop-stress  profile  in  figure  6  shows  the  steel  region  is  in 
compression,  resulting  from  the  winding  tension  of  composite  overwrap.  The  stress  is  not 
continuous  because  of  the  drastic  change  of  material  properties  (modulus)  from  the  composite  to 
the  steel  region.  It  also  illustrates  stress  relaxation  of  composite  overwrap  material. 
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4.  Relaxation  of  Thermal  Stresses 


The  time-dependent  thennal  viscoelastic  behavior  of  a  100-layer,  AS-4/3502  (graphite/PEEK) 
composite  cylinder  subjected  to  a  temperature  increase  AT  =  150°  C  is  examined.  Initial 
residual  stress  built  up  in  the  cylinder  due  to  the  AT.  The  composite  cylinder  has  an  inner  radius 
a  =  3.5  in,  an  outer  radius  b  =  4.1  in,  and  a  thickness  of  each  layer  h  =  6.0  x  10'3  in. 

Stacking  sequence  is  given  as  [0/3 0/60/90]25  from  inside  out,  with  the  0°  direction  coinciding 
with  the  axis  of  the  cylinder.  The  creep  property  of  an  AS-4/PEEK  graphite/epoxy  composite 
with  a  fiber  volume  fraction  of  0.67  was  measured  at  different  temperatures  by  Kim  and  Hartness 
(14).  The  study  shows  that  an  increase  of  compliance  with  time  due  to  creep  behaviors  of 
material  was  found  at  elevated  temperatures.  A  least-squares  curve  fitting  was  used  to  express 
the  transverse  and  shear  creep  from  the  original  AS-4/3502  data  in  power  law  forms  as  follows: 

S22 (t)  =  [ 1.7051  (i tf 1954  +  l]  S°22  (3?) 

and 

S66(t)  =  \l  1.3076  (i tf 2771  +  l]  Si  .  (38) 

Sr, 22  and  S°M  are  defined  in  the  previous  section.  The  compliance  in  the  fiber  direction, 

Sn  =  5.9  x  jo "V psi ,  S33(t)=  S 22(f),  and  Poisson’s  ratios,  v12  =  v, 3  =  0.3,v23  =  0.36 ,  are  assumed 
to  be  time-independent  for  this  example  simulation.  Thennal  expansion  coefficients  of  the 
composite  in  three  principal  directions  are  an  =  -  0.5  x  10'6  /°C  and 

a 2,  =  a33=  40.0  x  1  O'6  /°C  ,  where  the  negative  value  indicates  shrinkage  with  temperature 
increase. 

Figures  7  and  8  show  radial  displacement  and  radial  stress  profiles  across  the  thickness  of  the 
cylinder  at  three  instants — instantaneous  (initial  stress),  2  years,  and  infinite.  The  radial  traction 
and  displacement  satisfied  the  continuity  conditions  at  every  interface  of  layers  at  all  instants. 

The  radial  displacement,  w(t)  ,  will  reach  a  steady  state  over  a  long  period  of  time  (infinite  time) 

because  of  the  creep  behavior  of  composites.  In  fact,  the  radial  displacement  of  most  layers 
approaches  to  a  constant  value  except  at  the  innermost  and  outennost  portions  of  the  cylinder. 
The  free  traction  boundary  at  the  surface  of  cylinders  causes  the  gradients  in  the  radial 
displacements.  A  similar  phenomenon  is  also  observed  in  the  radial  stress  profile,  which 
approaches  to  a  constant  over  a  long  period  of  time.  This  long-term  creep  characteristic  reflects 
the  power  law  form,  equations  (21)  and  (22),  of  the  creep  compliance.  The  “saw”  shaped  radial 
stress  distribution  is  the  result  of  variation  of  fiber  orientations  through  the  thickness  of  the 
cylinder.  The  radial  stress  is  continuous  but  the  stress  gradient  is  not.  Accordingly,  the  stress 
profile  illustrates  the  “saw”  shape. 
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Figure  7.  Radial-displacement  profiles  in  a  cylinder  subjected  to  thermal 
loads. 


Figure  8.  Radial-stress  profdes  in  a  cylinder  subjected  to  thermal  loads. 


The  hoop  stress,  <x de(t) ,  through  the  thickness  of  the  cylinder,  is  illustrated  at  three  instants  in 
figure  9.  There  exist  two  distinct  values  (discontinuity)  of  oe9(t)  across  each  interface  of  two 
adjacent  layers  due  to  the  various  fiber  orientations  through  the  thickness.  The  hoop-stress 
profile  also  shows  a  trend  of  relaxation  over  a  period  of  time.  The  hoop-stresses  in  60  °C  and 
90  °C  layers  show  a  fairly  steep  gradient  across  the  cylinder  thickness  initially.  However,  the 
gradient  gradually  disappears  as  time  approaches  infinity. 
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Figure  9.  Hoop-stress  profiles  in  a  cylinder  subjected  to  thermal  loads. 


5.  Relaxation  of  Mechanical  Stresses 


In  the  following  study,  creep  and  stress  relaxation  of  a  composite  cylinder  subjected  to  internal 
pressure  is  investigated.  The  calculation  is  performed  using  a  0.6-in-thick  composite  cylinder 
with  a  lay-up  construction  of  [0/30/60/90]25-  Calculations  have  been  performed  using  similar 
basis  graphite/epoxy  composite  materials.  The  material  properties  are  the  same  as  described  in 
the  previous  section.  A  pressure  load  of  1000  psi  was  applied  at  the  inner  radius  of  the  cylinder. 
Some  selected  displacement  and  stress  profiles  from  the  model  prediction  are  illustrated  and 
discussed  in  the  following  sections. 

Figure  10  shows  the  radial  displacement  through  the  thickness  of  cylinder  at  three  instants — 
instantaneous  (initial  stress),  two  years,  and  infinite  time.  The  radial  displacement  profile  clearly 
illustrates  the  creep  behavior  of  cylinder.  The  inner  radius  of  the  cylinder  increases  over  a  period 
of  time  because  of  the  application  of  inner  pressure.  The  outer  radius  actually  shrinks  down 
because  of  creep  characteristics.  The  radial  strain  that  is  equal  to  the  gradient  of  radial 
displacement  increases  over  a  period  of  time. 

Figure  1 1  illustrates  the  relaxation  of  radial  stresses  through  the  thickness  at  three  instants.  The 
radial  stress  at  the  inner  radius  is  1000  psi,  equivalent  to  the  pressure  applied.  The  radial  stress  is 
zero  at  the  outer  surface  of  the  cylinder  since  it  is  traction-free.  The  curve  is  constructed  by 
connecting  the  stress  value  at  all  the  interfaces  of  layers.  The  “saw  shape”  of  curve  is  due  to  the 
variation  of  fiber  orientations  through  the  thickness.  Significant  relaxation  occurs  over  a  period 
of  time  as  shown  in  the  radial  stress  profiles. 
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Figure  10.  Radial-displacement  profiles  in  a  cylinder  subjected  to  internal 
pressure. 


Radius  (inch) 


Figure  11.  Radial-stress  profiles  in  a  cylinder  subjected  to  internal  pressure. 

The  hoop-stress  profiles  at  three  instants  are  illustrated  in  figure  12.  The  stress  is  not  continuous 
from  layer  to  layer  due  to  the  change  of  fiber  orientation.  The  gradient  of  the  stress  profile  is 
mainly  due  to  curvature  of  cylinder.  The  90  °C  layers  have  higher  stresses  because  they  are 
stiffer  in  the  circumference  direction  and,  thus,  carry  more  loads.  The  viscoelastic  effect  is  quite 
interesting  as  observed  in  the  hoop-stress  profiles  changes  over  a  period  of  time.  The  hoop  stress 
at  the  inner  radius  increases,  while  it  decreases  at  the  outer  radius.  The  integration  of  hoop  stress 
through  the  thickness  should  be  balanced  with  the  inner  pressure  applied  if  a  free  body  is  taken 
from  the  cylinder.  Since  the  stress  gradient  increases  over  a  period  of  time,  the  hoop-stress  will 
also  increase  at  the  inner  radius  of  cylinder. 
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Figure  12.  Hoop-stress  profiles  in  a  cylinder  subjected  to  internal  pressure. 


6.  Conclusions 


An  analysis  has  been  developed  for  viscoelastic  behavior  of  laminated  composite  cylinders  with 
ply-by-ply  variation  of  anisotropic  viscoelastic  properties,  which  cannot  be  studied  using  an 
isotropic  model.  Stress  relaxation  and  creep  are  properly  determined  in  a  cylinder  subjected  to  a 
thermal  and  mechanical  loads.  Emperimental  validation  was  conducted  using  a  high-tensioned, 
overwrapped  steel-ring  experimental.  Radial  stress  at  the  composite-steel  interface  was 
predicted  and  compared  with  the  measurement.  A  good  agreement  was  obtained.  The 
anisotropic  viscoelastic  behavior  of  the  composite  causes  interesting  characteristics  in  cylinders, 
which  are  critical  for  the  durability  of  the  structure.  Creep  and  stress  relaxation  could  exist  in  the 
fiber  direction  even  though  the  fiber  dominant  properties  are  elastic.  This  is  mainly  due  to  the 
contribution  of  the  Poisson's  ratio  effects  of  transverse  and  shear  properties.  Viscoelastic 
characteristics  are  critical  to  the  service  life  cycle  of  applications  such  as  pressure  vessels  and 
composite  rotors  designed  with  built-in  pre-stress  to  achieve  desired  mechanical  performance. 
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PRC  MAIL  CODE  R7000 
AUSTIN  TX  78712 

2  DEPUTY  AS  ST  SECY  FOR  R&T 
SARD  TT 

J  SINGLETON 
M  DONOHUE 
SAAL  TR  STE  9800 
2511  JEFFERSON  DAVIS  HWY 
ARLINGTON  VA  22202-391 1 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
RDRL SED 
E  SHAFFER 
BLDG  204  RM4C120 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

ABERDEEN  PROVING  GROUND 

62  DIR  USARL 
RDRL  SL 
R  COATES 
RDRL  VTM 
A  GHOSHAL 
DLE 
C  SHIAO 


RDRL  WM 
P  BAKER 
B  FORCH 
J  MCCAULEY 
RDRL  WML 
M  ZOLTOSKI 
RDRL  WML  D 
D  BEYER 
J  SCHMIDT 
RDRL  WML  F 
G  KATULKA 
D  LYON 
RDRL  WML  G 
D  DRYSDALE 
J  SOUTH 
RDRL  WML  H 
J  NEWILL 
B  SORENSEN 
RDRL  WMM 
J  BEATTY 
R  DOWDING 
P  SMITH 
J  ZABINSKI 
RDRL  WMS 
M  VANLANDINGHAM 
RDRL  WMM  A 
R  EMERSON 
L  HOLMES 
J  GARDNER 
M  GRIEP 
D  O’BRIEN 
T  PLAISTED 
J  SANDS 
E  WETZEL 
J  TZENG  (5  CPS) 

RDRL  WMM  B 
T  BOGETT1 
R  CARTER 
R  DOOLEY 
B  CHEESEMAN 
D  HOPKINS 
R  KARKKAINEN 
B  POWERS 
C  YEN 

RDRL  WMM  C 
R  JENSEN 
RDRL  WMM  D 
R  BRENNAN 
V  CHAMPAGNE 
ECHIN 
KCHO 
S  WALSH 
RDRL  WMM  E 
J  SINGH 
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NO.  OF 

COPIES  ORGANIZATION 


RDRL  WMM  F 
H  MAUPIN 
L  KECSKES 
RDRL  WMP 
S  SCHOENFELD 
RDRL  WMP  A 
J  POWELL 
RDRL  WMP  B 
C  HOPPEL 
S  SAT APATHY 
RDRL  WMP  D 
R  DONEY 
M  KEELE 
J  RUNYEON 
RDRL  WMP  E 
M  BURKINS 
P  SWOBODA 
RDRL  WMP  F 
A  FRYDMAN 
N  GNIAZDOWSKI 
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Intentionally  left  blank. 
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